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Abstract 

The surface evolution model based on geometric partial differential equation is used to numeri- 
cally study the kinetics of dewetting and dynamic morphologies for the localized pinhole defect in 
the surface of the ultrathin solid film with the strongly anisotropic surface energy. Depending on 
parameters such as the initial depth and width of the pinole, the strength of the attractive sub- 
strate potential and the strength of the surface energy anisotropy, the pinhole may either extend 
to the substrate and thus rupture the film, or evolve to the quasiequilibrium shape while the rest of 
the film surface undergoes phase separation into a hill-and- valley structure followed by coarsening. 
Overhanging (non-graph) morphologies are possible for deep, narrow (slit-like) pinholes. 

PACS numbers: 68.55.-a 
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I. INTRODUCTION 

Recent experiments with the sub-lOnm sihcon-on-insulator films at 800-900 °C 

demonstrate Si film dewetting caused by the long-range film-substrate interactions (which 
are also called wetting interactions). The mass transport in this system is by thermally- 
activated surface diffusion, and there is no stress at the film-substrate interface due to an 
absence of a lattice mismatch between the film and the Si02 substrate. Dewetting starts 
at randomly distributed pinhole defects in the Si planar surface. The pinholes may exist 
prior to the annealing, or they form shortly after the temperature is raised. Conditions 
favoring pinhole deepening over contraction, where the latter is caused by minimization of 
the surface area due to the mean surface energy (tension) , and the kinetics of the pinhole and 
its dynamic shape are not presently known. Besides, in contrast to liquids, the surface energy 
of solid surfaces is strongly anisotropic, leading to missing orientations in the dynamic or 
equilibrium surface shape i-B and facet., ..«ff-H.T.ese facto, are certain 
to affect the dynamics of the pinhole and moreover, through the nonlinear competition with 
the attractive wetting potential may lead to the emergence of an equilibria and thus to the 
suppression of the film dewetting and rupture. 

In Ref. IGj, following Ref. 19|, the partial differential equation (PDE)-based model is 
developed which allows to predict the wavelength of the fastest growing cosine-like perturba- 
tion of the film surface (also called the normal perturbation), assuming surface diffusion and 
the two- layer wetting potential 171]- 2l||. The model also enables computation of the dynam- 
ical (and probably faceted) morphologies. Such computations are performed for the normal 
perturbation and they demonstrate the stabilizing impact of the surface energy anisotropy 
on dewetting dynamics. The model can in principle support any reasonable form of the 



wetting potential 
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23| but the corresponding contribution to the governing PDE for the 



film thickness must be rederived. 

In this paper, using the model of Ref. [16|], the kinetics and morphologies are computed 
systematically for the localized surface defect from the full nonlinear PDE. As has been 
made clear above, real surface defects are necessarily localized. We compute for different 
widths and depths of the pinhole and, for all other model parameters fixed, observe very 
different dynamics and dewetting outcomes. We also relax the assumption made in Ref. 
[igI that the surface height above the substrate is described by a function h{x,t) (i.e., a 
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one-dimensional (ID) surface is nonoverhanging) and reformulate the model in terms o: 
parametric PDEs. This allows to compute, say, beyond the surface phase separation |8 
into orientations 0° and 90° for some surface energy anisotropics. Asymmetric morphologies 
and different kinetics may arise when the direction of the maximum surface energy is not the 
reference direction for the shape evolution (i.e., for instance, the ^-axis perpendicular to the 
substrate), which is often the case. Thus we incorporate such misorientation in the model. 



Note that parametric 
for instance Refs. 24 



brmulations of the geometric surface evolution laws are common, see 



25|,l26|. 



II. PROBLEM STATEMENT 

A 2D film with the free ID parametric surface T{x{u, t), z{u, t)) is assumed, where x and 
z are the Cartesian coordinates of a point on a surface, t is time and u is the parameter 
along the surface. The origin of a Cartesian reference frame is on the substrate, and along 
the substrate (x-direction, or the [10] crystalline direction) the film is assumed infinite. 
The 2;-axis is along the [01] crystalline direction, which is normal to the substrate. Marker 
particles are used to track the surface evolution (see Ref. 27|], for instance) thus x and z 
in fact represent the coordinates of a marker particle, which are governed by two coupled 
parabolic PDEs 
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Here the subscripts t and u denote differentiation, V is the normal velocity of the surface 



which incorporates the physics of the problem, and g = ds/du = \/ z\ is the metric 
function (where s is the arclength). 

The normal velocity of the surface is due to gradients of the surface chemical potential 
yU, which drive the mass flux of adatoms along the surface. In other words, redistribution of 
adatoms along the surface ch ang es its shape, which is equivalent to the surface moving in 

321 ]. The chemical potential is the sum of two contributions. 



the normal direction 



IsO, 3l|; 



[I = /i + where fi^'^^ is the regular contribution due to the surface mean curvature k, 



and /i^'") is the wetting chemical potential 

/i(-) = n i^,{9) - 75) — 
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20|: 



-z/i) 



cos 9, z > 0. 



(2) 



Here fl is the atomic volume, 6 is the angle that the unit surface normal makes with the [01] 
crystalline direction, 75 = const, is the surface energy of the substrate in the absence of the 
film, i is the characteristic wetting length, and 7p(6') is the primary part of the anisotropic 
surface energy of the film, i.e. for typical four-fold anisotropy 



7(^) = 70(1 + e, cos4{e + P)) + -K^ ^ 7p(^) + i^^^' 



(3) 



In Eq. ([3]) 70 is the mean value of the surface energy, determines the degree of anisotropy, 
j3 is the misorientation angle and 5 is the small non-negative regularization parameter having 
units of energy. The 5-term in Eq. ([3]) makes the evolution equations ([1]) mathematically 



m 

, [ll| . (The anisotropy is weak when < 



well-posed for strong anisotropy [33| - [35| 

< 1/15 and strong when e^>l/15. 5 = in the former case.) The surface energy has a 
maximum at 4(^ + /?) = 0, i.e. at ^ = —d. For^ = 0, this direction is the 2;-axis. 



The curvature contribution to /i is 
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-2;/^))K + 75exp(-2;/£)/t \^^^s 
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(4) 



where the subscript s denotes differentiation with respect to the arclength. If the wetting 
potential is zero [zjt ^ 00), this reduces to familiar strongly anisotropic form js, [o], 11 1. 
Finally, 



and 



(5) 



(6) 



where D is the adatoms diffusivity, v the adatoms surface density, k the Boltzmann constant 
and T the absolute temperature. The only differences of this formulation from the formu- 



lation in Ref. 



16|, except for accounting for the surface energy misorientation in Eq. ([3]) 



and the parametric representation, are that in the wetting (exponential) contributions 
are accounted for in full (that is, the approximation in the form of averaging across the film 
thickness is not employed) and the regularization term is not included in ^^"^K (/i^"") does 
not contain the surface stiffness 7 + 7e6i and thus it does not make the PDE ill-posed for 
strong anisotropy. Besides, the regularization contribution to /i^"') is vanishingly small for 
large surface slopes due to its proportionality to h'"^ [16]). 

To nondimensionalize the problem, the thickness of the planar undisturbed film, h^, is 
chosen as the length scale, and hi/ D as the time scale. Also, let r = £//io- The dimensionless 



problem is comprised of Eq. ([5]) and Eqs. ([T]) (where differentiations are with respect to the 
dimensionless variables), and where 

V = B {ftif + /ii:)) , (7a) 
/^^'^^ = (jp + ^) (1 - i-'M) «^ + rexp {-z/r)K - A (y + ^^s.}j , (7b) 

= (7,(e) - r) ^^Piz^M cos e, (7c) 
7p(0) = l + e^cos4(^ + /5). (7d) 

In Eqs. ([7]) i? = Q'^U'^o/^kThl), T = 'js/lo and A = 6/{'johl). For the computational 
method, using the relation between s and u, 

- = -- , (8) 

the problem is written entirely in terms of the independent variables u and t (not shown). 

In the simulations reported below, the following values of the physical parameters are 

used: D = 1.5 x 10~^ cm^/s, = 2 x 10"^^ cm^, 70 = 10^ erg/cm^, 75 = 5 x 10^ erg/cm^, 

u = 10^^ cm"^, kT = 1.12 x 10^^^ erg, Hq = 10^^ cm, and 5 = 5 x 10"^^ erg. These values 

translate into B = 3.57 x 10^^, F = 0.5 and A = 5 x 10^'^ (or zero). In this paper we 

consider strong anisotropy, > 1/15 and r = 0.02, 0.1. 

The initial condition in all runs is the Gaussian surface 

2" 



z{x, 0) = 1 — dexp 



X 



w 



< X < 10 (9) 



where < d < 1 and w are the depth and the "width" of the pinhole at t = 0, respectively 
(see Figure 4). The length of the computational domain equals to ten times the unperturbed 
film thickness, and the defect is positioned at the center of the domain. We use values 
d = 0.5, 0.9, which correspond to shallow and deep pinhole at t = 0, respectively, and 
w = 0.15 (narrow pinhole), w = 1 (intermediate pinhole) and w = 2 (wide pinhole). 

The method of lines is used for the computation with the periodic boundary conditions 
at X = 0, 10. Eqs. ([I]) are discretized by second-order finite differences on a spatially- 
uniform grid in u. The integration in time of the resulting coupled system of the ordinary 
differential equations is done using the implicit Runge-Kutta method. Initially u = x, but 
periodically (usually after every few tens of the time steps) the surface is reparametrized 
so that u becomes the arclength, and the positions of the marker particles are recomputed 



accordingly. This prevents marker particles from coming too close or too far apart in the 
course of surface evolution. 

III. RESULTS 
A. Kinetics 

Figures 1 and 2 show the log-normal plots of the pinhole depth vs time, for d — 0.9 and 
d — 0.5, respectively. Zm is the height of the surface at the tip of the pinhole. 

Wide and intermediate deep pinholes dewet but the depth of the narrow deep pinhole 
decreases until it reaches quasiequilibrium at 2; = 0.75 (Figure 1). Quasiequilibrium means 
that Zm (or, cquivalently, the depth) changes very slow or not at all, while the rest of the 
shape may change relatively fast. Correspondingly, we will call the surface shape at the time 
when the quasiequilibrium depth is attained, the quasiequilibrium shape. (Again, this shape 
is changing, but the height of its minimum point (the tip) is not changing, or is changing 
very slow). Also notice in Figures 1(a) and 1(b) that the growth rate at the rupture time is 
finite for r = 0.1 but infinite (or extremely large) for r = 0.02, and the time to rupture is 
somewhat less for r = 0.1. However, from Figure 1(c), the time to reach quasiequilibrium is 
about ten times larger for r = 0.1 than for r = 0.02. 

As can be seen in Figure 2, only the wide shallow pinhole dewets, and only when r = 0.1. 
In all other cases of w and r (except w — 1, r = 0.1, shown by the solid hne in Figure 2(b)) 
the quasiequilibrium is achieved. For w — 1, r = 0.1, the depth is initially a non-monotonic 
function of time, but after the transient phase it monotonically and slowly decreases without 
reaching the quasiequihbrium (we computed for i < 2 x 10^). Also, one can see that for 
r = 0.02 the depth at quasiequilibrium decreases as w decreases, and it takes less time to 
reach quasiequilibrium as w decreases. 

In Figure 3 the dewetting kinetics is compared for several misorientations and strengths 
of the anisotropy. Shown is the case of the deep, wide pinhole and r = 0.1. The time to 
rupture increases insignificantly with the decrease of anisotropy or with the increase of the 
misorientation angle. The faster dewetting for stronger anisotropy here can be attributed 
to the initial faster shape changes due to larger gradients of i.e. before the surface 
orientation falls into an unstable (spinodal) range and faceting steps in, and to the proximity 
of the pinhole tip to the substrate. That the misorientation slows the kinetics is well-known. 
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for instance Liu & Metiu |8|] in their important study of faceting call the similar situation an 
"off-critical quench" , and find that at sufficiently large misorientations the "crystal surface 
will not phase-separate spontaneously, but will have to overcome a finite free-energy barrier" . 



B. Morphologies 



Figure 4 shows the initial and the final surface shapes of the initially deep pinhole for the 
three values of w. The kinetics of the corresponding dynamical shapes is shown in Figures 
l(a)-(c) by solid lines and has been discussed above. Dewetting of the wide and intermediate 
pinholes proceeds through the extension of the tip of the pinhole until it reaches the substrate 
at 57°. The quasiequilibrium shape for the narrow pinhole is very similar to the one shown 
in Figure 5(c). The latter shape is discussed in more details below. While evolving from 
the initial slit-like shape to the quasiequilibrium shape, the surface of the narrow pinhole 
overhangs, until it slowly returns to the non-overhanging shape later (Figure 5). In the time 
interval where the overhanging takes place, the surface slope is large and non-analytic. 

Surface shapes for the initially shallow pinhole are shown in Figure 6. Characteristic 
of these shapes is emergence of the hill-and- valley structure 5|], which becomes possible 
even when the film dewets (case of the wide pinhole in Figure 6(a)) due to larger time 
required to reach the substrate in this case. The angle at rupture is 90°. Figure 6(b) shows 
coarsening of the structure for the intermediate width case. As has been pointed out above 
at the discussion of Figure 2(b), in this case there is no quasiequilibrium (at least until 
t > 2 X 10'^), and the tip recedes towards the unperturbed height h = 1. Note that the 
apparent recession rate is slower than the overall coarsening rate, t = 7240 is the time when 
the two pyramidal structures appear on the film surface. The slope of their walls is shown in 
the inset and it is almost constant in each of the four characteristic intervals of x, with each 
interval corresponding to a facet. For 7240 < t < 2 x 10'^ the walls (facets) become more 
straight and the graph of z^ix) becomes nearly constant in each of the four characteristic 
intervals. Figure 6(c) shows the quasiequilibrium shape for the narrow pinhole. Formation 
of the hill-and- valley structure followed by coarsening continues for t > 4, as is evidenced by 
the surface slope shown in the inset at the left, but the pinhole depth at x = 5 is constant. 
(In fact, the difference of depths at t = 197 and at t = 4 is 0.007.) Note that formation of the 
hill-and-valley structure, its coarseningand slope selection have been the subject of many 



papers, see for instance Refs. 



B-Q, H, G 



20j, (l7|, |36|-|4^. (Refs. 



13, 



m 



20| discuss 



impacts of wetting interactions). Since in this paper we are interested in characterizing 
dewetting and rupture, we do not further pursue that direction. 

Finally, Figure 7 demonstrates impacts of misorientation (/? = 10°) on morphology. As 
expected, asymmetrical shapes emerge for /5 7^ 0°. While the wide, shallow pinhole dewets, 
only its right sidewall undergoes the phase separation into a hill-and-valley structure. Ki- 
netics is very similar to the case [3 = 0, see Figure 2(a). 

IV. DISCUSSION 

In this paper a fully nonlinear model is used to compute the complex scenarios of dewet- 
ting/equilibration for a localized pinhole defect in the surface of strongly anisotropic thin 
solid film, assuming the two-layer wetting potential and zero lattice mismatch with the 
substrate. 

The computed dewetting kinetics can be explained at large using the magnitude of the 
dewetting factor. In Ref. 3] it was shown that the dominant dewetting terms in the 
mass-conservation evolution PDE are proportional to exp(— z/r)/r^ (for r < 1). This is 
plotted in Figure 8. One can see that for r = 0.1, 0.02 the maximum of this factor occurs 
for small values of the film thickness. This explains why all but one initial conditions in 
Figure 1 (deep pinhole) lead to pinhole depth increase until the film dewets, while all but 
one initial conditions in Figure 2 lead to decreasing pinhole depth. The dewetting factor can 
also explain the small (large) difference in characteristic time scales for r = 0.1 and r = 0.02 
in Figure 1(a) (2(a)). Indeed, for the case of Figure 1(a) the factor is 36.8 (r = 0.1) and 
16.8 (r = 0.02) vs, respectively, 0.7 and 3 x 10"*^ for the case of Figure 2(a). Evolution of 
the film also depends strongly on the film shape and whether 6 is in the unstable range. 
When it is, as in the cases shown in Figures 4-7 then, the faceting instability is energetically 
more favorable than dewetting and most often the competition of the two processes causes 
the unusual hill-and-valley structure where shallow pinhole remains despite the structure 
coarsening (which takes place separately at both shoulders of the pinhole). Note that the 
faceting instability seems to be always initiated where the surface changes from horizontal 
to a sidewall (regions A and B in Figure 4). 

We also point out that in our computations the pinhole tip is always nonfaceted at rupture 
(see Figures 6(a) and 7(a)). In fact, dewetting ceases if the facet spreads to the tip. Also, if 
the pinhole is not too deep (Figure 4), then the contact angle with the substrate at rupture 
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is 90°. In reality, after the contact this value decreases in order to minimize the total energy 
of the surface-substrate system 
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FIGURES CAPTIONS 



Figure 1. Kinetics (rate) data for deep pinhole {d — 0.9). e-y = 1/12, (5 — 0. Line slope 
equals the rate of the tip evolution. Sohd line: r = 0.1. Dash line: r = 0.02. (a) Wide 
pinhole {w — 2). (b) Intermediate pinhole {w — 1). (c) Narrow pinhole {w — 0.15). 

Figure 2. Same as Figure 1, but for shallow pinhole {d — 0.5). 

Figure 3. Kinetics data for deep, wide pinhole, r = 0.1. Sohd line: e-y = 1/12, P — 0. 
Dash line: e-y = 1/12, P — 10°. Dash-dot hne: = 1/8, P — 0. Dash-dash-dot-dot line: 

= 1/14, p = Q. 

Figure 4. Surface morphologies for the deep pinhole, r = 0.1, = 1/12, /3 = 0. The 
depth (d) and the width (w) of the pinhole at t = are defined. Solid, dash, dash-dot line: 
w — 2, 1, 0.15, respectively. Inset: magnified view of dewetting region for w = 2, 1. 

Figure 5. Magnified view of the transient surface morphology and surface slope (inset) for 
the deep, narrow pinhole from Figure 4. Dash line: surface at t = 0. 

Figure 6. Surface morphologies for the shallow pinhole, r = 0.1, — 1/12, P — 0. (a): 
w — 2, (b): w — 1, (c): w — 0.15. Dash line: surface at i = 0. 

Figure 7. Surface morphology for the shallow, wide pinhole, r = 0.1, e-y = 1/12, P — 10°. 
Dash hne: surface at t = 0. 

Figure 8. Contour plot of the dewetting factor exp {—z/r)/r'^. 
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